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Abstract 

We propose a method, based on a Generalized Heisenberg Algebra 
(GHA), to reproduce the anharmonic spectrum of diatomic molecules. 
The theoretical spectrum generated by GHA allows us to fit the ex- 
perimental data and to obtain the dissociation energy for the carbon 
monoxide molecule. Our outcomes are more accurate than the stan- 
dard models used to study molecular vibrations, namely the Morse 
and the q-oscillator models and comparable to the perturbed Morse 
model proposed by Huffaker pQ, for the first experimental levels. The 
dissociation energy obtained here is more accurate than all previous 
models. 

Keywords: molecular vibrational spectrum; Heisenberg algebra; quantum- 
group. 
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1 Introduction 



Vibrational spectroscopy has a fundamental role in molecular physics and its 
applications extend to other fields like astronomy [2|, biology |3] and earth 
and environmental sciences jlj. Vibrational molecular analysis provides im- 
portant informations on the structure of the molecules. Recently, the de- 
velopment of the powerful experimental techniques which allow the study of 
highly excited vibrational states jSHH]- I n Born-Oppenheimer approximation, 
electronic, rotational and vibrational quantum states [3j can be considered 
separately and molecular vibrations can are approximately described by the 
harmonic oscillator [7|. However, this model fails to describe the highest 
molecular vibrational modes and consequently to provide the correct disso- 
ciation energies. One of the first attempts to improve these approximations 
is due to Morse jH] who introduced a simple potential which allows to solve 
exactly the Schrodinger Equation (SE) and provides a reasonably approxima- 
tion to the spectrum of diatomic molecules, including an upper bound which 
lacks in the harmonic oscillator potential. Later on, many other approxima- 
tions based on the use of Morse-like, Kratzer-like or modified versions of these 
or other potential functions have been introduced in order to obtain better 
fittings with experimental data EH HU E2] and for evaluating dissociation 
energies jT31 E] . 

Although ab initio methods have shown to be able to provide good ap- 
proximations for the spectra of many molecules, some authors have pointed 
out some of their limitations. Angelova et al. [15 , noted that the Morse 
potential contains only quadratic corrections and one needs to use the em- 
pirical Dunham expansion to fit the highest vibrational energies. Iachello 
and Levine [16J, remarked that the solution of the SE is very difficult in 
the case of two- and three- dimensions and, at this point, algebraic methods 
could bring some advantages. Problems related to the dissociation energy es- 
timation and also a algebraic energy method to calculate it, were described 
by Sun et al. [T7]. Nevertheless, besides the facts above mentioned, there 
are still some molecules whose potentials deviate strongly from Morse-like 
potentials [TS] . 

The algebraic approach to the study of molecular vibrational and rota- 
tional spectra was pioneered by F. Iachello in jTH] (for a recent review see 
|20j). The method has been proved to be useful for calculating accurately 
highly excited vibrational levels and for describing the whole spectrum of 
complex molecules, while maintaining its simplicity in cases where the use of 
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ab initio methods are not feasible in practice |51 I23|. 

Deformed Heisenberg Algebras (DHA or q-oscillators) have attracted con- 
siderable attention since they have been proposed [2U due mainly to 
their potential physical applications (see ref. [27] for a review). They 
are non-trivial generalizations of the Heisenberg Algebras characterized by 
the deformation parameter q. It is interesting to notice that for some values 
of the q- deformation parameter, namely for < q < 1, the energy spectrum 
associated with the Deformed Heisenberg Algebras presents an upper bound, 
as in fact occurs in systems like composite particles [2E], as for instance, di- 
atomic molecules ^SJ |2H HUB EI]- Within this formalism g-oscillators are 
used as a model of anharmonic oscillators, and the value of q is chosen to 
reproduce experimental molecular vibrational spectra. When compared with 
other methods, the DHA approach requires only a few parameters, allows 
analytical expressions and fits well the first experimental levels. However, as 
it will become clear later, it does not provide always the correct dissociation 
energy because it may fail to fit the highest vibrational levels. Recently, 
many authors have discussed the physical interpretation of this parameter 
[221 1221 1211 EH] but this analysis is out of the scope of this work. 

Recently, a Generalized Heisenberg Algebra (GHA), where the commu- 
tation relations among the operators depend on a characteristic function of 
the generalized number operator has been proposed - If the characteristic 
function is linear, with slope q 2 , the algebra corresponds to the g-oscillator 
algebra. In the present work we introduce and implement a nonlinear GHA 
which is able to describe typical features of the vibrational molecular spec- 
trum of the carbon monoxide molecule (CO). The spectrum generated by 
this nonlinear GHA (nl-GHA) allows us to fit the 20 first vibrational tran- 
sitions and to obtain the correct dissociation energy, providing better global 
results than the methods above mentioned. 

The paper is organized as follows. In section El we discuss the solution 
of the Schroedinger Equation (SE) with Morse potential . In section 3 we 
introduce the GHA method. In section 4 we apply our GHA model to the 
CO molecule. Finally, section 5 is devoted to discussions and conclusions. 
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2 Morse Potential for Diatomic Molecules: 
Vibrational Levels 

Let us now turn our attention to one of the most familiar quantum descrip- 
tions of diatomic molecules. For low energy levels, the molecular potential 
can be approximated by the harmonic potential, whereas for high energy 
this approximation breaks down. In fact, the higher vibrational transitions 
exhibit a certain degree of anharmonicity. Thus, to take into account anhar- 
monic terms and rotational effects, nonlinear terms must be introduced in 
the energy expression jSj: 

Euj = J2 = Y ^(^+l) l [AJ+i)r, (1) 

l,m ^ ' 

where v and J are the vibrational and the rotational quantum numbers, 
respectively. The above expression was obtained by Dunhan, using the WKB 
method [Tj. 

An analytical solution for the SE, which corresponds to the second order 
truncation of the equation ((Tj) with J = 0, was obtained by Morse [S] using 
the so-called Morse potential (see Figure UJ): 

V(r) = De~ 2a{r ~ To) - 2De- a(r - ro) , (2) 

where a is a characteristic constant of each molecule, r is the nuclear distance, 
tq is the equilibrium position and D is the dissociation energy. The solution 
of the SE for this potential gives the following expression for the vibrational 
energy 

e n = -D + hu[(u + i) - X e{y + ~) 2 ], (y = 0, l, (3) 

where Xe is a positive constant called the anharmonicity coefficient, u) is the 
fundamental frequency, v is the vibrational quantum number and k — — > 1. 
The anharmonicity coefficient is related to the dissociation energy by, 

Xe = ^ (4) 

The energy levels obtained agree reasonably well with experimental data 
for many molecules [3J. The maximum allowed value of n in equation (JHJ) is 
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such that the SE solution remains finite. One can also show that this value 
is given by 

11 

Hence, there is maximum energy value e Umax beyond which the system is not 
bounded. 

In order to improve the Morse results, several perturbative methods, 
based on the Morse or other analytical or numerical potential, have been 
proposed. Huffaker [1], used a perturbed Morse model to calculate the first 
terms of the Dunhan expansion (jl} for the CO molecule. He used a potential 
given by following expression, 

V(r) = V e {{\ - e^f + 6 4 (1 - e- a ") 4 (6) 
+b 5 (l - e- ac >) 5 + ...)]. 

Using 28 spectral lines of the CO molecule [37|, he calculated 8 parameters of 
the potential curve given by equation (j7J). This model is good when compared 
with the model obtained via RKR method by Mantz et al. [37] ■ If one ignores 
the rotational degree of freedom (1 = 0), the model has 7 parameters. 



3 The Generalized Heisenberg Algebra 

In j^H] Curado e Rego-Monteiro proposed an algebra called GHA which is 
generated by three operators Jo, A and A\ satisfying the following relations: 

J ^ = ^ /(J ), (7) 
A Jo = /(J ) A, (8) 
[A A ] ] = /(Jo) -Jo, (9) 

where t is the Hermitian conjugate, (At)t = A, Jq = J and /(Jo) is an 
analytical function of Jo, called the characteristic function. Assuming the 
existence of a vacuum state represented by |0) (A\0) = 0), it can be shown 
that 

J Q \m-l) = f m -\eo)\m-l),m= 1,2.., (10) 
At|m-1) = M m _i|m), (11) 
A\m) = Af TO _i|m-l), (12) 
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where M^ l _ 1 = / m (e ) — e , e is the lowest J eigenvalue and / m (e ) is the 
m-th iteration of the function /. This algebra describes a class of Heisenberg- 
like algebras of quantum systems, having energy eigenvalues given by 

£n = /(Cn-l), (13) 

where e n and e n _i are two successive eigenvalues and are related with the 
energy for the system [38J. Therefore, the Jo eigenvalues can be obtained 
iteratively and they can be upper bounded or not, depending on the chrac- 
teristic function, the values of the function parameters and the initial value 
€q. For each kind of function, the values of the parameters and the initial 
values determine the existence or not of the fixed points, e* = /(e*) and 
their stability. As a consequence, different spectra are obtained for different 
functions and the eigenvalues behavior can be analyzed using dynamical sys- 
tems techniques. For the linear case, /(Jo) = Q Jo + s, the equations (JZHHI) 
become 

[J ,Al] q = sA\ (14) 
[Jo,A] rl = —A, (15) 

[A,A*] = (q-l)J + s, (16) 

where [a, 6] = ab — qba is the g-commutation of the operators a and b. 
Equations (jl4H16j) describe the one-parameter GHA studied in jHSj and can 
be mapped into the g-oscillator algebra. If q = 1, the standard Heisenberg 
algebra is recovered. 

A graphical analysis of the functions /(e) = qe + s and y = e is shown 
in Figure El The intersection between the two lines is identified as the fixed 
point of the recurrence equation e n = ge n _i + s. The most interesting cases 
are obtained for < q < 1 and e < y^, since under these conditions the 
energy spectrum has an upper bound e* = jj-, as observed in the spectra of 
bounded systems, as for instance diatomic molecules 12Z] ■ 

4 GHA and Vibrational Molecular Spectrum 
4.1 Linear case - (g-oscillator) 

We will now study the spectrum generated by the linear GHA mentioned 
above, which correspond to the g-oscillator algebra. To study the vibra- 
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tional spectrum of diatomic molecules via GHA formalism, we use a result 
developed in [2211301 • Let us start with the general Hamiltonian 

H = Huj(ciAA^ + c 2 A^A + c 3 ), (17) 

where A and A^ obey the relations (JZEl) and c\, c 2 and C3 are real numbers. 
Choosing c\ — c 2 — 1 and C3 = 0, and using the equations (jlOH12j) we get [SH] 

H = hw(f(J ) + J -2e ). (18) 

Therefore, for the linear case /(Jo) = Q^o + s (q < 1) we have: 

H = Hu((q + l)J + s-2e ). (19) 

Applying the Hamiltonian (fT9*|) on the eigenstate |z/) of Jo (with h — 1) we 
obtain iJ|^) = E u \u), where the energy eigenvalues are given by 

E u = u[(q+l)f v (e ) + s-2eo], (20) 
and f u {eo) = q v t§ + S ^~T' After some algebra we obtain 

E v = to (M q - L 9 <f +1/2 ) , (21) 

with L q = i+f (e - 1^) = ^7§(eo - e*) and M q = 2(^ - e ) = 2(e* - e ). 

The maximum energy [v — > 00) for this system is given by E^ = uM q = 
2tu(e* — e ), and the dissociation energy is then 

D q = E oo -E = Lu(l + q)(e*-e ). (22) 

The q— parameter can be interpreted as being related to the anharmonicity 
constant, Xe = — | r , where q = e T ffifj . 

We will now apply the previous results to the carbon monoxide molecule 
spectrum and compare it with both experimental data and the spectra ob- 
tained by using the Morse, the perturbed Morse and the g-oscillator models. 
We used the data provided by the HITRAN database jH] and chose the 
lines in which the molecules are in the electronic fundamental state. Also, 
we selected the lines whose rotational quantum numbers are zero. Thus, 
there remained the spectral lines which represent the first 21 pure vibra- 
tional transitions of the CO molecule. 
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In order to fit the parameters (q and eo) to the experimental data, note 
that the logarithm of the difference between two successive levels of the 
equation (|2ip. as a function of u, is a straight line with slope Inq. Taking the 
logarithm of the difference between two successive data points, we obtain the 
value of q by evaluating the slope of the curve after a linear regression. The 
fitted parameters are shown in table (0). The linear GHA spectrum agrees 
well with the experimental data (see table (J3j) and figure (J3J)) mainly in the 
case of the 20 first vibrational transitions. However, the D value calculated 
by equation (|22j) deviate strongly from the experimental value as can be 
seen in table (0). This means that, although the g-oscillator relative errors 
(AE = Ethe °£~ Eex p- ) are smaller than the relatives error obtained by Morse 

model, for most of the 20 first vibrational transition (figure (J5J)), the approach 
based on the g-oscillator is not able to fit higher frequency vibrational levels. 



4.2 Nonlinear case 

As we have seen above, the linear GHA (i.e. the q-oscillator) , is a good 
model only for the first levels. In order to fit the energy levels and to obtain 
the correct dissociation energy we were led to use a nonlinear functional f(x) 
in the GHA approach (J7EJ). We found that, by using the functional 

f{x) =p x 4 + q x + 1, (23) 

it is possible to obtain a good fit with the experimental data and the correct 
dissociation energy. Due its nonlinearity, the spectrum generated by this 
nonlinear GHA does not have an analytical closed expression like equation 
(121)) . In this case the spectrum can only be calculated numerically. 
For the sake of simplicity, we start from the Hamiltonian 

H = hw{A ] A + e ) = hwJ , (24) 

which is derived from Hamiltonian (JT7j) with c\ — 0, c-i — 1 and C3 = eo- 
Replacing /(Jo) given by equation (j2l?|) in relations (JTEl) and applying H on 
the eigenstates of Jo (with h = 1) we obtain 

E u = wr(e ), (25) 

where f v = e v is given by: 

e n+1 = pe A n + qe n + l. (26) 
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We use the values generated by equation (|25|l with f u (e ) given by equation 
(|26|) in order to fit the q, p and eo parameters with the experimental data. 
The dissociation energy is given by 

D GHA = E oo -E = w{e* -e ), (27) 

where e* is the stable fixed point of the recurrence equation (J26j) . 

The fitted parameters are shown in table (jlj). The dissociation energy 
are shown in table (J2J). In table (JSJ) we compare the energy levels obtained 
with the different models referred in this work with experimental data. The 
relative errors are shown in the figure (JSJ). We can see that the nonlin- 
ear GHA with the functional (|23|) provides better fittings with experimental 
data than both q-oscillator and Morse models and it is comparable to the 
Huffaker model [I]. Furthermore, nonlinear GHA provides a more accurate 
dissociation energy when compared with that obtained in all other methods. 
This results show that our model is a good method of obtaining the higher 
anharmonic energy levels of the CO molecule. 

As the p parameter is very small (table ([TJ), if eo is also small, the first it- 
erations of the equation (J2T)|) are dominated by the linear term. The nonlinear 
term becomes relevant as the number of iterations increase. The difference 
between two successive energy levels decreases up to zero as the function 
()26|) is iterated. Because the p parameter is negative, the difference between 
two successive energy levels tends to zero. Thus, the nonlinear term flats the 
energy curve faster than in the pure linear and the perturbed Morse cases 
but slightly slower than the Morse Model, as we can see in the figure (@J). 

We would like to stress that besides these accurate fittings, the nonlinear 
GHA gives us an extremely simple way to estimate higher transitions for the 
CO molecule, just by iterating equation (|2*U|) up to the required level. 

5 Conclusion and Perspectives 

In this work we have proposed a method based on GHA to reproduce the 
vibrational molecular spectrum of diatomic molecules. For the CO molecule, 
we have shown that the previous method based on the (/-oscillator algebra 
(linear GHA) reproduce the vibrational molecular spectrum for the 20 first 
vibrational transitions, but it does not provide the correct dissociation en- 
ergy, i.e., this model fails when describing higher vibrational levels. Using 
a nonlinear functional (4 — th order) in the GHA we were able to fit the 
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experimental energy levels and also to calculate the dissociation energy for 
the CO molecule with good accuracy For this molecule the nonlinear GHA 
has provided better global results than the usual Morse, perturbed Morse 
models and (/-oscillator model. This work shows that GHA can be used as a 
phenomenological simple tool to study composite particles. 

GHA produces different spectra for different characteristic functions. The 
GHA spectrum could be related to deformations in harmonic potential. In- 
deed, this algebra could reproduce a variety of spectra once we are able to find 
the appropriate characteristic function. Consequently, this algebra could be 
used to study other molecules, simply choosing the appropriate characteristic 
function. In an even more general fashion, de Souza et al. [I2| have con- 
structed a more general structure that depends on two functionals. This new 
algebraic structure can reproduce a even greater variety of spectra allowing 
us to work with two quantum numbers. 

The results obtained in this work motivate us to further explore the re- 
lation between the GHA and molecular systems. This investigations may 
include: 1) the study of the relation between the functional parameters and 
the molecular parameters; 2) the investigation of the potential curve under- 
ling the GHA and 3) the account of electronic and rotational states. 
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model 


parameters 


values 


Morse 


Xe 


0.0062 


(/-oscillator 
(linear GHA) 


q 


0.98646 


eo 


36.98 


nonlinear GHA 


eo 


0.9235 


Q 


0.9872 


V 


-1.43 x 10~ y 


perturbed Morse 
(Huffaker) 


o 


77.21317 


T 


83769.28cm- 1 


64 


0.036067 


h 


0.017505 


h 


0.014945 


h 


0.010770 


bs 


0.008142 



Table 1: Values of the parameters used in each model. 
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model 


Dissociation Energy (cm r ) 


Experimental 


89591.35 


Morse 


86426.44 


(/-oscillator 


158970.48 


Perturbed Morse 1 


95394.23 


Perturbed Morse 2 


95394.23 


nonlinear GHA 


89987.76 



Table 2: Values of the dissociation energy calculated for each model. 



transition 
line 


energy (cm r ) 


exp. 


TV T 

Morse 


1 /"I T T A 

nl-GHA 


1 • /"I T T A 

linear GHA 


i n 1 t\ t '-< 

perturbed Morse d 


1 — > 


2143.24 


2143.30 


2144.56 


2152.46 


2143.27 


2 -> 1 


2116.76 


2116.39 


2117.10 


2123.32 


2116.79 


3^2 


2090.34 


2089.48 


2089.99 


2094.57 


2090.37 


4^3 


2064.00 


2062.57 


2063.19 


2066.21 


2064.02 


5^4 


2037.72 


2035.66 


2036.68 


2038.23 


2037.73 


6^5 


2011.51 


2008.75 


2010.44 


2010.63 


2011.51 


7^6 


1985.38 


1981.84 


1984.43 


1983.41 


1985.35 


8^7 


1959.32 


1954.93 


1958.62 


1956.55 


1959.26 


9^8 


1933.33 


1928.01 


1932.97 


1930.06 


1933.24 


10 -> 9 


1907.43 


1901.10 


1907.43 


1903.93 


1907.29 


11 -> 10 


1881.61 


1874.19 


1881.99 


1878.15 


1881.40 


12 -> 11 


1855.85 


1847.28 


1856.58 


1852.72 


1855.59 


13 -> 12 


1830.19 


1820.37 


1831.18 


1827.63 


1829.84 


t 14 -> 13 


1804.61 


1793.46 


1805.74 


1802.89 


1804.16 


15 -> 14 


1779.11 


1766.55 


1780.23 


1778.48 


1778.55 


16 -> 15 


1753.69 


1739.64 


1754.61 


1754.39 


1753.00 


17-> 16 


1728.36 


1712.73 


1728.85 


1730.64 


1727.53 


18 -> 17 


1703.12 


1685.82 


1702.90 


1707.21 


1702.13 


19 -> 18 


1677.96 


1658.91 


1676.75 


1684.09 


1676.79 


20 -> 19 


1652.88 


1632.00 


1650.37 


1661.29 


1651.53 



Table 3: Vibrational spectrum of the CO molecule. Comparison of the exper- 
imental data with the Morse, ^-oscillator, perturbed Morse and GHA model. 
Experimental data from HITRAN database. 
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Figure 1: Morse (full lines) and harmonic potential (dashed lines) energy 
levels. The horizontal axis is the relative nuclear position r. Note that the 
two potential have a minimum value \h r = r$ and that the Morse Potential 
has a dissociation energy D. 



f(e) 




Figure 2: Graphical analysis of a linear characteristic function, /(e) for < 
q < 1 with e < y£^. Note that, as the function is iterated, the energy value 
approaches the stable point e* = yz^ (the upper bound). The dissociation 
energy of the system is proportional to (e* — e ). 




Figure 3: Comparison between the relative errors of the Morse (o), q- 
oscillator (o), nonlinear GHA (□) models and perturbed Morse (*). We 
can see that the nonlinear GHA errori^ire smaller (in almost all data points) 
than the Morse and g-oscillator models. 




20 40 60 80 100 120 

vibrational quantum number 



Figure 4: Normalized and shifted energy spectrum generated by the Morse 
(*), g-oscillator (□) and nonlinear GHA (•) models, perturbe Morse (o) and 
experimental data (+). The dashed l?ri^ represents the experimental dissoci- 
ation energy value. We can see that the dissociation energy of the nonlinear 
GHA lies on the experimental dissociation energy (dashed curve) while the 
g-oscillator and Morse spectrum are respectively above and below dissocia- 
tion energy line. Inset figure shows in detail the Morse and the nonlinear 
GHA energy curves close to dissociation energy line. Note that the Morse 
curve flats before reaching the exnerimental dissociation enerev while the 



